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_ Well determined radial velocities and abundances are essential for analyzing the properties of the Globular Cluster 

, system of the Milky Way. However more than 50% of these clusters have no spectroscopic measure of their metallicity. 

I In this context, this work provides new radial velocities and abundances for twenty Milky Way globular clusters which 

lack or have poorly known values for these quantities. The radial velocities and abundances are derived from spectra 
(N- obtained at the Call triplet using the FORS2 imager and spectrograph at the VLT, calibrated with spectra of red 

^ , giants in a number of clusters with well determined abundances. For about half of the clusters in our sample we 

\l ■ present significant revisions of the existing velocities or abundances, or both. We also confirm the existence of a sizable 

■ abundance spread in the globular cluster M54, which lies at the center of the Sagittarius dwarf galaxy. In addition 
evidence is provided for the existence of a small intrinsic internal abundance spread ((T[Fe/H]i„t « 0.11-0.14 dex, 

, similar to that of M54) in the luminous distant globular cluster NGC 5824. This cluster thus joins the small number of 

■ Galactic globular clusters known to possess internal metallicity ([Fe/H]) spreads. 

^ ■ Key words. Stars: abundances - Stars: kinematics and dynamics - Stars: Population II - Galaxy: globular clusters - 

■ Galaxy: globular clusters: individual: NGC 5824 - Galaxy: stellar content 
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1. Introduction where relative ages were calculated for 64 clusters. With 

respect to metallicities, Carretta et al. (|2009p assembled a 

Galactic globular clusters (GGC) represent one of the fun- ^^^^ jjj ^^j^^^ ^33 dusters: the objects were put 

daniental systems that allow a reconstruction of the early ^ ^-^^^^ ^^^j^^ metallicities were computed based 

evolution of the Milky Way: if their ages, metallicities and j^^j^g^ published in four different studies. Twenty five 

kinematics were known with sufficient precision, then cor- ^j^^^^^g ^^^^ ^^1^^^ ^^^^ ^^^-^ resolution spectroscopic 

relations of their kinematics and chemical abundances with ^^^^^ Carretta & Gratton (1997, CG97) and Kraft and 

time would shed light on the dynamical and chemical eyo- ^^^^^ (pgggi) - which have 15 objects in common, and the 

lution of the protogalactic halo and bulge. To reach this ^^^^ ^^^^ 2inn & West ZW84) and Rutledge 

goal, large and homogeneous data samples are needed, and ^^gg^^^^ j^g^^ r^^^ -^^^^ 2W84 is based on 

indeed considerable progress has been seen m recent years .^^g^^^^g^ i^gf^f^ narrow-band imaging, so the largest ho- 

The largest collectioii of color-magnitude diagrams (CMD) j^ogeneous spectroscopic study of individual globular clus- 

obtamed with a single telescope (HST) and umform data ^^^^^ ^^-^ ^j^^^ j^g^ However their sample repre- 

reductions is that of Marm-Franch et al. (EDMl MF09), ^^^^^ 44% ^^jects in the Harris (2010) cata- 
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logue, which is a shortcoming for many investigations of 
the kind illustrated above. For example 17 clusters of the 
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MF09 study (~ 25% of the total) do not have an entry in 
R97, and their metaUicity was taken essentially from ZW84 
(after a transformation to the CG97 scale). Another exam- 
ple is offered by Forbes & Bridges (|2010l) who compiled 
ages and metallicities for 93 GGCs, and based on their age- 
metallicity relation (AMR) they estimated that ~ 25% of 
the clusters were accreted from a few dwarf galaxies. Most 
of their table entries come from MF09, but for 29 clusters 
ages and metallicities had to be extracted from inhomoge- 
neous sources. R97 based their work on spectra obtained 
with the 2.5m Dupont telescope at Las Campanas, so the 
absence of homogeneous data affects mostly outer-halo (or 
heavily extincted) distant clusters, which is where relics of 
accretion are most likely to be found. 

To remedy this situation, we commenced a project to 
collect medium-resolution spectra for clusters not included 
in the R97 sample. Spectra in the z-band region of giant 
stars were obtained with F0RS2 at the ESO/VLT obser- 
vatory to measure metallicities based on absorption lines 
of the Can triplet (CaT). We employ the "reduced equiv- 
alent width" method first introduced by Olszewski et al. 
p^grj) and Armandroff & Da Costa (fTMTl hereafter AD91), 
which is the same as used in R97. One caveat with this 
method is that a-elements (which include Ca) are enhanced 
in GGCs compared to the solar value, with [a/Fe]'^ -1-0.3 
(e.g. Carney I1996p . so it is important to identify clusters 
with anomalous enhancements, especially when choosing 
the calibrators of the index-[Fe/H] relations. For this rea- 
son, we also collected F-band spectra to explore a metaUic- 
ity ranking based on the strength of Fe and Mg lines. For 
this spectral region, we are in the process of deriving effec- 
tive temperature Tes, gravity log 5 and metallicity [Fe/H] 
from a full spectrum fitting of sample stars using a library 
of 1900 stars observed at high resolution with the ELODIE 
spectrograph, as described in Katz et al. (|201ip . The results 
will be presented elsewhere (Dias et al., in preparation). 

In a first paper based on our data (Da Costa et al. l2009p 
we discovered a metallicity spread in the giant stars of M22 
almost in parallel with the high resolution study of Marino 
et al. (|2009p . Here we publish the full catalog of CaT-based 
metallicities and radial velocities, comprising 20 clusters. 
As already discussed, the main point of our project is the 
homogeneity of the data and analysis, through which we 
add ~ 30% more objects to the R97 sampl^B- Furthermore, 
our calibration is based on a very large set of GC templates 
accurately chosen to have well established abundances from 
high-resolution studies. 

A few more points can also be gathered from Table [1] 
where for each observed cluster we list the "best" source 
of metallicity that we could find in the literature. For two 
clusters (NGC 6139 and NGC 6569) the metaUicity is still 
based on the integrated-light Q^g index, and in the case 
of NGC 5824 it is based on the integrated-light version 
of the CaT method (Armandroff & Zinn il988)) . Because 
individual stars are not measured, for these clusters pos- 
sible internal metallicity spreads cannot be discovered, as 
the case of NGC 5824 clearly shows (see below). In nine 
cases high- resolution studies exist, but almost all by dif- 
ferent authors, causing concerns on the relative abundance 
ranking. In addition high-resolution studies require large 
observational and data analysis investments, limiting the 



^ We had planned observations for 49 clusters, but poor 
weather conditions did not allow completion of the program. 
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Figure 1. A selection of typical fully reduced spectra is 
shown, one for each of the 28 clusters observed. For each 
spectrum the cluster and star identification are given above 
the continuum. 

number of stars that can be measured and thus the sta- 
tistical significance of metallicity dispersions, if they are 
found. Finally six clusters have [Fe/H] estimated via the 
equivalent width of Fe lines measured on medium-resolution 
spectra, and calibrated on high-resolution spectra of stan- 
dard clusters. Again different methods and calibrations are 
applied by different authors, raising concerns on the self- 
consistency of the metallicity scales. To summarize, there 
are only two clusters, where the same CaT method em- 
ployed here was used (NGC 2808 and Pyxis) to determine 
abundances. Hence we expect that the (reduced) EWs pub- 
lished here will be the basis for all future abundance ranking 
of the clusters listed in Table \T\ 

2. Observations and reductions 

2.1. Selection of targets 

The core target list was assembled as follows. We started 
from the catalog of GGCs published in Harris (jl996l ver- 
sion Feb. 2003, hereafter H96) and considered all objects 
visible during ESO P77, i.e. having RA between 12 h and 
h and DEC southern of -1-20° (which makes them observ- 
able from Paranal for at least 2hr at airmass X < 1.5). We 
also required to reach at least My = — 1 along the RGB 
(see AD91), and to have 1/ = 20 as the faintest RGB star 
to limit exposure times to ~ 1 hr. Further, to exclude the 
closest clusters, a requirement (m — M)v > 16.75 was also 
imposed. Finally we used the compilations of Pritzl et al. 
(2005) and Gratton et al. (|2004p to remove those clusters 
that had CaT measurements and to add calibration clus- 
ters. A few objects that we deemed interesting for a number 
of reasons were then added back to the list:Terzan 7 of the 
Sagittarius dwarf spheroidal galaxy, NGC7006 as 'second 
parameter' object, and NGC6325, NGC6356, NGC6440, 
NGC6441, and HPl as bulge members. The final list even- 
tually had 57 entries, including the 8 calibrators. 
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Table 1. Observations log 



Cluster UT tcxp (sec) notes N [Fe/H]' 



May 28, 2006 



NGC3201 


23:06 


2x7.5 


c 


17 




M4 


03:02 


2x1.9 


c 


16 




MIO 


03:16 


2x20.9 


c 


19 




NGC6397 


05:31 


2x2.6 


c 


17 




NGC6528 


08:44 


2x61.2 


c 


19 




NGC6553 


09:04 


2x31.1 


c 


18 




M71 


09:23 


2x6.1 


c 


14 




M15 


09:39 


2x18 


c 


18 




Ruploe 


23:22 


2x349 


2nd 


15 


h/BWZ97 


NGC2808 


00:14 


2x16.8 




19 


CaT/R97 


NGC5824 


00:35 


2x248.9 




19 


i/AZ88 


Lynga? 


01:26 


2x200.2 




19 


Fe/TF95 


NGCd139 


03:34 


IxIbAA 




19 


i/Q39 


Terzan3 


03:55 


2x230.5 




19 


Fo/C99 




U4:i5 


2x214.1 


b 




Fe/M95a 


NGC6356 


04:36 


2x69.4 


b 


19 


Fc/M95a 


NGC6380 


04:59 


2x615.6 




19 


Fo/C99 


NGC6440 


05:48 


2x295.4 


b 


19 


h/OVR08 


NGC6441 


06:37 


2x96 


b 


19 


h/OVR08 


NGC6558 


07:03 


2x60.7 


b 


19 


h/B+07 


Pal7 


07:25 


2x97.1 




18 


Fe/C99 


NGC6569 


07:51 


2x88.3 




19 


i/Q39 


M22 


08:14 


2x13.3 


m. 1 


19 


h/BW92/- 


M22 


08:28 


2x13.3 


m. 2 


19 


-/LBC91/M+09 


M22 


09:55 


2x13.3 


m. 3 


18 








May 29, 2006 






Pyxis 


23:23 


2x1258.9 


2nd 


18 


CaT/PKMOO 


HPl 


05:31 


2x488.1 


b, X 


39 


h/B+06 


M54 


08:15 


2x278.8 




19 


h/BWG99/- 












-/B+08/Carl0 


Terzan? 


08:53 


2x334.9 


sag, X 


25 


h/T+04/S+05 


NGC7006 


09:20 


2x680 


2nd, X 


25 


h/K+98 



Notes. N is the number of stars observed in each cluster. 
c=calibration, X=done with mask (MXU) to increase slitlets 
density, b=bulge, sag=Sagittarius dSph, 2nd=second parame- 
ter. The metallicity source is coded as method/references, where 
the reference acronym can be found in the bibliography list. The 
method can be (h)igh resolution, (i)ntegrated light, medium- 
resolution EW of (Fe) lines, or (CaT) method. 

2.2. Observations 

Our data were obtained with F0RS2 (Appenzeller et al. 
I1998p . working at the Cassegrain focus of VLT/UTl-Antu. 
The instrument offers two ways of realizing multi-object 
spectroscopy, either with 19 movable slitlets (MOS mode), 
or by fabricating masks to be inserted in a mask exchange 
unit (MXU mode). To define the masks and MOS slitlets 
configurations, nine hours of pre-imaging in the V and / 
bands were used to observe the 49 main targets in ser- 
vice mode, while pre-imaging for the calibrators was avail- 
able from previous VLT runs. The 6'.8 x 6'.8 field-of-view 
was centered on the cluster for the sparser or more dis- 
tant systems, but was offset from the centers for the nearer 
and richer clusters. For each cluster the resulting color- 
magnitude diagram was used to select targets on the RGB. 



Afterward, spectra were collected in visitor mode, during a 
two-night run at the end of May 2006. The choice of MXU 
or MOS mode was based on the stellar density. The CaT 
spectral region was covered with the 1028z+29 grism and 
the OG590+32 order-blocking filter, yielding a maximum 
spectral coverage of -7700A-9500A at a scale of 0.85A per 
(binned) pixel. In all cases the slit width was defined at 1". 
For the MOS mode the number of slits is sometimes smaller 
than 19, because of the targets' distribution on the sky (see 
Table [T]). As the table shows, in three cases masks were 
used to cover higher density clusters, with up to 39 slitlets 
per mask. The target star magnitudes typically cover '^3 
mag along the RGB, and two exposures were obtained to 
allow removal of cosmic-rays. Because of the limited time 
and bad weather conditions, spectra could be collected for 
20 clusters only (plus the eight calibrators). The observing 
log is reported in Table [TJ 23 clusters could be observed 
in the first night, while clouds in the second night limited 
observations to 5 clusters only. 

2.3. Extraction of spectra 

All spectra were extracted using the F0RS2 pipeline ver- 
sion 1.2 (Izzo & Larsen 120081) . which splits the processing 
into two steps. First, using daytime calibration frames, slit 
positions are found, and the wavelength calibration and dis- 
tortion map for each slit are created. Then the products of 
the first step are used to reduce the target spectra. The 
pipeline logs give a mean residual scatter around the wave- 
length calibration of 0.24 pixels (i.e. 0.2 A or 7kms^^ for 
a dispersion of 0.82 A px~^), a mean spectral resolution of 
« 2440, and a mean fwhm of the arc lines of 3.51 ± 0.07 A 
(130 kms~^). Given this resolution we expect individual 
radial velocities to have uncertainties of 10-15 kms^^, and 
thus the mean velocity for a cluster will have an uncer- 
tainty 3 to 4 times smaller depending on the number of 
members. This is consistent with the uncertainties given in 
Table [21 We note though that a comparison of our mean 
cluster velocities with well determined values in the litera- 
ture (see Sect. 13. ip suggests that the true uncertainties in 
our mean velocities are actually somewhat larger, probably 
as a result of less easily quantified systematic effects such 
as mask-centering errors, etc. 

The software corrects for bias and fiatfield, and com- 
putes a local sky background for each slit, to be subtracted 
from the object spectra. The Horne (1986) optimal extrac- 
tion is applied. The spectra are normalized by exposure 
time, and the wavelength solution is aligned to a refer- 
ence set of >20 sky lines by applying an offset (0.09 px 
on average). The two exposures were average-combined af- 
ter pipeline processing. The S/N ratios for the final spectra 
varied from ^110 for the brightest stars to ~25 for the 
faintest in a typical exposure. To illustrate the quality of 
our data, Fig. [1] presents one spectrum for each of the ob- 
served clusters. 



3. Radial Velocities 

Heliocentric radial velocities of the single stars were com- 
puted with the rvidlines task in iraf's rv package using 
all three of the Can triplet lines. Afterward, the cluster 
velocity was computed by averaging the velocities of the 
cluster members. In doing so, stars with significantly dis- 
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Table 2. Radial Velocities 



o 
o 



Cluster 


N ( 


vr) ±err 


(vr) ±err 








kms 


-1 


kms 


1 


kms^^ 


HIO 


a) metal-poor standard clusters 








iVi4 


1 K 


/ Z 


6 


70.7 


0.2 


1.3 


iViiU 


1 Q 
lo 


Q« 
oO 


2 


75.2 


0.7 


10.8 


In 'orODoy ( 


ID 


1 9 


3 


18.8 


0.1 


-6.8 


M71 


11 


— 9^ 


4 


-22.8 


0.2 


-2.2 


M15 


-LO 


-111 


5 - 


107.2 


0.2 


-4.0 




17 


486 


4 


494.0 


0.2 


-8.0 


n 1 TTi t^'i~ CI 1 - 


JJUUi JJiL 


^rarr 


clusters 






Pyxis 


8 


42 


4 




1.9 


7.6 


NGC2808 


17 


91 


4 


101.6 


0.7 


-10.6 


RuplOe 


9 


-39 


4 


-44.0 


o.u 


o.u 


NGC'iS'?^ 

I'l vjr ij<J 


17 


-11 


2 


-27.5 


1.5 


16.5 


Lynga 7 


9 


22 


3 


8.0 


o.u 


o.u 


NGC6139 


15 


34 


4 


6.7 


6.0 


27.3 


Terzan 3 


13 - 


-131 


4 - 


136.3 


0.7 


5.3 


NGC6325 


10 


28 


2 


29.8 


i.o 


— I.o 


xlr i 


8 


41 


4 


45.8 


0.7 


-4.8 




5 - 


-198 


4 - 


197.2 


1.6 


—0.8 


Pol 7 


14 


164 


3 


155.7 


1.3 


8.3 




7 


-47 


3 


-28.1 


5.6 


-18.9 


M22 


51 - 


-150 


2 - 


146.3 


0.2 


-3.7 


M54 


19 


137 


2 


141.9 


0.5 


-4.3 


NGC7006 


20 - 


-379 


1 - 


384.1 


0.4 


5.1 


c) metal— 


rich standard clusters 








NGC6528 


4 


205 


2 


206.6 


1.4 


-1.6 


NGC6553 


18 


-9 


4 


-3.2 


1.5 


-5.8 


d) metal— 


rich prof 


;ram 


clusters 








NGC6356 


12 


67 


4 


27.0 


4.0 


39.0 


NGC6380 


9 


7 


3 


-3.6 


2.5 


10.6 


NGC6440 


8 


-76 


4 


-76.6 


2.7 


1.6 


NGC6441 


7 


18 


4 


16.5 


1.0 


1.5 


Ter7 


14 


154 


3 


166.0 


4.0 


-12.0 



Notes. AwR is in the sense FORS2-H10. HPl: The dif- 
ference listed is with respect to the heliocentric velocity given 
by Barbuy et al (2006), which we give in lieu of the Harris cat- 
alog value. It appears the catalog lists the uncorrected observed 
value. 



crepant line strengths were excluded even if they had ve- 
locities compatible with cluster membership. The velocity 
errors were computed as standard errors of the mean as 
determined from the standard deviation of the individual 
velocities. Our results are summarized in Table [2] where N 
is the number of cluster members used in calculating the 
mean velocity. 

In detail our assessment of cluster membership for each 
star observed was based on two assumptions: that the range 
in the radial velocities of member stars was small, and that 
the dispersion in the measured equivalent widths about the 
fitted line was comparable to the measurement errors. The 
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Figure 2. Membership selection for stars observed in 
NGC6569. The upper panel shows the heliocentric radial 
velocity plotted against V — Vhb- The dotted line is the 
cluster velocity given by the Harris catalogue, the dashed 
line is the mean velocity for the 7 stars ultimately selected 
as cluster members. The stars unambiguously identified as 
radial velocity non-members are plotted as x-signs. The 
lower panel shows the sum of the strengths of the A8542 
and A8662 lines plotted against V — Vhb- Here, 4 stars, 
plotted as open circles, have line strengths incompatible 
with the remaining 7, plotted as filled circles, which have 
a small dispersion about a line of slope -0.627 A mag~^ 
(see Sect. U]). These 4 stars are classified as line-strength 
non-members. 



latter assumption is equivalent to assuming the intrinsic 
abundance dispersion in a cluster is small, an assumption 
that despite recent discoveries (e.g.. Sec. I5.1.3|) remains 
generally valid: most globular clusters are mono-metallic. 
Nevertheless, careful consideration was given to the likely 
membership in every case where an observed star had a ve- 
locity consistent with that of the cluster but a line strength 
different from that expected for its V — Vhb value, given the 
line strengths of other candidate members. For the more lu- 
minous stars in the more metal-rich clusters a discrepant 
weaker line strength was often the result of a depressed 
pseudo-continuum caused by the ~8440 A TiO bandheads 
(see Olzsewski et al. 1991). Such occurrences were read- 
ily identifiable and the stars showing TiO were included in 
the cluster velocity calculation, if the velocity was consis- 
tent with other members, but not in the determination of 
the cluster line strengths. The final selection of members 
was also required to define a sensible sequence in the color- 
magnitude diagram derived from the cluster pre-imaging, 
though given the process by which the stars to observe were 
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initially selected, this provides only a consistency check. We 
also note that the spectra of every star observed was visu- 
ally inspected so that any effects from the occasional poorly 
removed cosmic ray or inadequate sky-subtraction could be 
allowed for. 

We illustrate our membership selection process in Fig. [2] 
for the cluster NGC6569 which is typical of the sample of 
program clusters. The upper panel shows the measured ra- 
dial velocities for 18 of the 19 stars observed in this cluster. 
The spectrum of the remaining star had an instrumental 
defect in the vicinity of the A8662 line of the Ca II triplet 
and so was not used in the subsequent analysis. The figure 
shows that there are 7 stars whose velocities differ signifi- 
cantly from the others; such stars are unambiguously non- 
members. The lower panel shows the sum of the strengths of 
the A8542 and A8662 lines of the Can triplet plotted against 
V — Vhb- Here four of the stars with velocities compatible 
with cluster membership have line strengths that are sig- 
nificantly larger than those for the remaining seven, whose 
line strengths are consistent with measurement error driven 
scatter about a single (line strength, V — Vhb) relation. 
Consequently, we consider these 4 stars as line strength 
non-members - the alternative that there is a second pop- 
ulation in the cluster with metallicity ~0.8 dex higher is 
much less plausible. 

We also verified that our velocity results are not affected 
by the choice to measure the line centers of the three Ca ii 
triplet lines rather than using cross-correlation techniques. 
Specifically, we reanalyzed the velocities for the standard 
cluster NGC 3201 and the program clusters NGC 5824 and 
HPl by applying the f xcor task within iraf. The spec- 
tral region correlated was 8400-8770 A which is largely free 
of atmospheric absorption contamination and which con- 
tains a number of weak stellar lines. For NGC 3201 and 
NGC 5824 the template employed was a cluster member 
while for HPl the template used was MIO star 1_3150. 
The mean difference between the velocities from the cross- 
correlation process and from the triplet line center mea- 
surements was 0.8 km s^^, 0.8kms~^ and l.Okms^^ for 
NGC 3201, NGC 5824 and HPl, respectively The standard 
deviations and the number of stars used were 5.6kms~^ 
(16), 1.8kms-i (17) and 4.1kms-i (26), respectively. For 
7 HPl member stars the mean difference and standard devi- 
ation were l.Gkms^^ and 3.9kms~^, respectively. We con- 
clude that our adopted velocities and their uncertainties are 
not affected by the choice of measurement technique. 

3.1. Comparison with Harris catalog values 

To verify the zero-point of our velocities, we first com- 
pared our velocities with those from H96 (further updated 
in December 2010, hereafter HIO) using only the clus- 
ters with a listed uncertainty of less than 2kms~^. There 
are 20 such clusters, and for these the mean difference 
is VpoRS — Vmo — 0.10 km s~^ with a standard devia- 
tion of 7 kms~^. This indicates the true uncertainties in 
our mean velocities are likely of order 5-6kms~^ rather 
than 3-4kms~^, the mean of the errors listed in Table [2] 
Nevertheless, the effectively zero mean offset indicates ex- 
cellent consistency between our determinations and existing 
well determined ones. The velocity differences Ffors — V'hio 
for these clusters are shown as open circles in Fig. |31 

The largest difference among these clusters is for 
NGC 5824 where we find a velocity of -11 ± 2 kms^^ 




-400 -200 200 

Heliocentric Velocity (FORS) 



Figure 3. The velocity differences between our values and 
those of HIO are plotted here versus our velocities. The 
open circles are for clusters where Harris lists an error of 
less than 2kms~^ and filled stars are for clusters where the 
HIO listed error exceeds 2kms~^. The vertical error bars 
are given by the quadratic sum of our and HIO errors. 



from 17 members while the Harris catalog lists a veloc- 
ity of —27.5 ± 1.5 kms^^. The catalog entry is essen- 
tially that of Dubath, Meylan & Mayor p997p who give 
Vf — —26.0 ±1.6 km s~^ from an integrated spectrum of the 
cluster centre. Earlier less precise values, e.g. —30 kms~^ 
(Armandroff & Zimii988J, -28kms-i (Zi nn fc W est rTOM)) 
and —38 kms^^ (Hesser, Shawl & Mever [T986I hereafter 
HSM86) are not inconsistent with the Dubath et al. (|1997p 
value. This suggests that our NGC5824 velocities might be 
all systematically ~15-20 kms~^ too high. Such offset rep- 
resents 1/8 to 1/6 of the FWHM resolution element, and 
could be due to an overall mis-centering of the MOS mask. 
This interpretation is supported by a comparison of our 
velocities with those determined from high resolution spec- 
troscopy of three NGC5824 stars by Villanova and Geisler 
(private communication). The mean velocity found from 
the high dispersion spectra is — 30.8±2.4 kms^'^ while for 
the same three stars our mean velocity is — 9±3 kms~^, 
supporting the idea of a systematic offset in the F0RS2 
velocities. We emphasize though that the existence of this 
possible velocity offset does not affect in any way the classi- 
fication of the 17 F0RS2 stars as NGC 5824 members. We 
intend to obtain additional spectra of stars in this cluster, 
which will either confirm or lead to a revision of our value. 

For the remaining eight clusters in our sample, the un- 
certainty listed in the HIO catalogue exceeds 2 kms~^. The 
difference between our new determinations and the cata- 
logue values are plotted as filled stars in Fig. |3l For five 
of the clusters the previous estimates are not inconsistent 
with our newer and better established values, but for the 
remaining three (NGC 6139, NGC 6569 and NGC 6356) 
the discrepancy with the previous estimates is substantial. 
For NGC 6139 and NGC 6569 the only previous deter- 
minations come from relatively low resolution spectra. In 
particular, for NGC 6139 the previous determinations are 
8 ± 7 kms^i from Webbink pMTj) and 4 ± 12 kms^^ from 
HSM86. Our value of 34 ± 4 kms~^ is undoubtedly prefer- 
able. For NGC 6569 our determination (—47 ±3 kms~^) is 
in reasonable accord with that given by Zinn & West (|1984p . 
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—36 ± 14 kms^^, but disagrees with that, —26 ± 6 kms^^ 
given by HSM86. Again our new determination is to be 
preferred. 

The largest discrepancy is for NGC 6356 where we find 
a velocity of 66 ± 4 kms~^ from 13 member stars while 
the catalog of HIO lists a velocity of 27 ± 4 kms~^. Our 
value is supported by the value given by Minniti (|1995ap 
of 56 ± 5 kms^^. It is possible that in generating the cata- 
log entry the (labeled uncertain) HSM86 velocity has been 
over weighted and the Minniti (1995) value under weighted. 
Regardless of the origin of the discrepancy it is likely that 
our value is preferable. 

Finally Terzan 7 and HPl deserve further comments. 
The velocity quoted by HIO for the first cluster is that of 
DA95 (4 stars), which is 12kms~^ higher than our value. 
However Sbordone et al. (|2005p have published velocities 
for 5 stars that are in better agreement with our result. 
Averaging the numbers from their Table 1, the velocity is 
159 ± lkms~^, with a difference of only +5kms~^ com- 
pared to us. In the case of HPl, it appears that the velocity 
quoted in HIO is the geocentric value given in Barbuy et al. 
(|2006p . while the heliocentric value agrees with our result 
within the errors. Indeed that velocity is the one given in 
the HIO column of Table H 

4. Metallicities 

Metallicities were computed by following the well estab- 
lished method originally proposed by AD91. Initially, for 
each star in all clusters the equivalent widths of the two 
strongest CaT lines (A8542, A8662) are measured, and their 
sum T,W is computed. Then using only calibration clus- 
ters, 'EW values are plotted vs. V — Vhb, and the slope 
of the T,W vs. V — Vhb relation is found (see Fig. H]). To a 
good approximation the slope a is independent from metal- 
licity within the range spanned by GGCs, so it can be 
used to define so called reduced equivalent widths W' — 
1^8542 + M^8662 — a {V ~ Vhb) for all stars in each cluster. 
In this way the gravity dependence of EWs is empirically 
removed, and an average value (W) can be computeqj- 
Finally the [Fe/H] vs. (W) relation for the calibration clus- 
ters was fit with a low-order polynomial to define the cal- 
ibration relation (see Fig. [S]) . At this point W' values can 
be computed for each star in the rest of the clusters, and 
they can be converted into [Fe/H] values by means of the 
relation found above. 

To obtain a more reliable abundance calibration, the 
standard cluster observations from Gullieuszik et al. (|20091 
G09) were also analyzed. These spectra were obtained with 
the same instrument setup and were reduced using the 
same methods employed here, the only difference being 
the slit width (O'.'S for G09 and 1" here). The additional 
calibration clusters from G09 are NGC 4590, NGC 4372, 
NGC 5927, NGC 6397 (no stars in common with our data 
set), NGC 6528 (no stars in common with our dataset), 
NGC 6752, M5 and NGC 6171. The total of calibration 
clusters is then 14 (see Table [5|). 

It should be noted that the CaT method works on the 
assumption that Vhb depends mainly on [Fe/H], and that 
other stellar parameters play a secondary role. While this is 
true for the relative ages of old stellar systems like globular 

^ (W) also coincides with the intercept of the linear fit to the 
{T,W,V - Vhb) points 



cluster^, cluster-to-cluster differences in helium abundance 
might instead be significant, and might introduce additional 
scatter in the metallicities yielded by the method. In gen- 
eral GGCs share a common He abundance (Buzzoni et al. 
nM31 ZoccaU et al. [MM Cassisi et al. [MJ5|) . but things 
might be different for bulge clusters. Recently Nataf et al. 
(12011 1 explained the difference between the luminosity of 
the RGB bump of the Galactic bulge, and that predicted by 
the luminosity-metallicity relation of Galactic globular clus- 
ters, by postulating that bulge stars have an He enhance- 
ment Ay = 0.06 (see also Renzini IT^M)) . To see the effect 
of such enhancement on the metallicity computed via the 
CaT method, one would need to take stellar models com- 
puted with different helium abundances, and then compute 
EWs of Ca lines and V-band luminosities using atmospheres 
with gravities and effective temperatures that are changed 
accordingly. Indeed we plan to carry out these tests in a 
forthcoming paper, while here we can check what is the ef- 
fect on the luminosity of the HB of the He enhancement 
quoted above. According to Renzini (|1977p the bolometric 
luminosity of the HB varies as AMhb = —4.7 x AY at 
a fixed metallicity and age. The bulge HB, and presum- 
ably that of its clusters, might therefore be brighter by 
AVhb ^ 0.28 mag. Such shift in luminosity would shift 
(W) by 0.28 X 0.627 = 0.18 A, where 0.627 is the slope 
of the EW, V - Vhb relation (Fig. g]). The change in (W) 
causes a change in [Fe/H] that depends on the metallicity 
itself (Fig. El), being ~ 0.05 dex for metal-poor clusters, and 
up to ^ 0.2 dex for metal-rich ones. However bulge clusters 
in Fig. [5] have very small [Fe/H] differences with respect 
to the literature, and in line with those of the rest of the 
clusters. An exception is HPl, but the difference with the 
literature is not systematic: our [Fe/H] is smaller than that 
of HIO and larger than that of C09/Appendixl by the same 
^ 0.3 dex, which is more than what an He enhancement 
would predict. The conclusion is that the cluster-to-cluster 
scatter in He content should not affect the results of this 
work. 

In the following sections the analysis is described in 
more detail, starting with a description of how photomet- 
ric catalogs and equivalent widths were obtained, and then 
moving to the [Fe/H]-(W) calibration relation. 

4.1. Photometry 

The pre-imaging data (which included both short and 
long exposures) were used to create photometric catalogs. 
Stetson's daophot/allstar package fStetson ll987lll994() 
was used to carry out the Point-Spread-Function photome- 
try. Instrumental magnitudes were calibrated by using color 
terms and zero points provided by ESO as part of their rou- 
tine quality controQ. We estimate that the zero-points on 
the standard V, / system are uncertain at the 0.05-0.10 
mag level because the preimaging observations were not 
all taken in photometric conditions. This is not a big con- 
cern however, because the method uses relative luminosi- 
ties. Photometry of G09 clusters was taken from previous 
work, with the exception of NGC6528 (see below). 



^ The effect of age on the CaT method is discussed in Da 
Costa & Hatzidimitriou (1998), Cole et al. ([MM)l . Pont et al. 
p004]) . and Carrera et al. 1 2007)1 . 

* see http: / /www. eso.org/ obser ving/dfo/quaUty /F0RS2/ qc / 1 
/ photcoeff/ photcoetfs fors2.html 
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Figure 4. Plot of Ca II line strength (^3542 + ^^8662) against magnitude difference from the horizontal branch {V — Vhb) 
for the standard clusters. In order of increasing ^8542 + ^8662 values at ^ — Vhb = -1-5, the solid lines are for clusters 
M15 (individual stars plotted as open 6-pt star symbols), NGC 4590 (open triangles), NGC 4372 (filled triangles), NGC 
6397 (open 5-pt stars), MIO (open circles), NGC 6752 (filled 5-pt stars), NGC 3201 (filled diamonds), M4 (x-symbols), 
M5 (plus symbols), NGC 6171 (filled squares), M71 (open diamonds), NGC 5927 (open squares), NGC 6528 (6-pt star 
symbols), and NGC 6553 (filled circles). The data for each cluster has been fit with a line of slope — 0.627A/mag for 
V — Vhb < 0.2. Vertical bars on each point show the measurement uncertainty in the line strengths. 



4.2. Measuring equivalent widths 



Equivalent widths were measured for the two strongest CaT 
lines in the co-added spectra by fitting a model profile over 
the line central bandpasses as defined by AD91, and then 
by computing the encompassed area. In AD91 Gaussian 
functions were used, which provide an excellent fit to CaT 
line profiles for clusters with metallicities up to 47 Tuc 
([Fe/H] ~ —0.7 dex). On the other hand it was realized 
by Rutledge et al. (|1997b|) and Cole et al. (,2004. C04) that 
CaT lines of stars in more metal-rich clusters have pro- 
files that are not Gaussian because of the strong damping 
wings, and they proposed to fit Moffat functions or sums of 
Gaussian and Lorentzian functions. In this work we tried 
to have an independent view on this problem, and we used 
calibration clusters to test both Gaussian and Gaussian 
plus Lorentzian (G+L) functions. Gaussian fits were per- 
formed as described in AD91, and we followed C04 for the 
G+L fits. These adopted a common line centre A™ and 
the best-fit parameters were determined using a Levenberg- 
Marquardt least-squares algorithm (see Markwardt I2009|) . 
To estimate the equivalent width measurement errors, we 
also measured the Ca 11 line strengths independently on the 
individual spectra of each star. The calibration clusters of 
G09 were already fit with G+L functions, so only Gaussian 
fits were performed in this case, to obtain the correspond- 
ing TiW values. The resulting 'EW obtained with the two 
methods are compared in Appendix [Aj where we show that 
a linear transformation exists between the two sets of mea- 
surements. We therefore defined G09 as our reference YiW 



scale, and transformed into this system the widths com- 
puted with the AD91 method. 

Having verified that there is a one-to-one correspon- 
dence between widths measured with any of the two 
methods, we measured metal-rich program clusters with 
G+L fits, and for metal-poor program clusters SVF were 
computed with Gaussian fits and transformed into the 
G09 scale. The metal-poor group contains clusters Pyxis, 
Rupl06, NGC5824, NGC6139, Ter 3, NGC6325, HPl, 
NGC6558, NGC6569, M22, M54 and NGC7006, while the 
metal-rich group contains clusters NGC6356, NGC6380, 
NGC6440, NGC6441, and Ter 7. In addition NGC2808, 
Lynga 7 and Pal 7 were measured with both methods to 
further confirm the reliability of the transformation defined 
in Appendix [XI 

Coordinates, radial velocities, V — Vhb values, and 
equivalent widths for the cluster member stars are pub- 
lished in tables to be found in the electronic version of the 
paper (see Tables [3] and 2] for an example of the layout). 
These are the fundamental measurements of this work, 
which allow a different metallicity calibration to be applied 
in the future, should it become available. 

4.3. Reduced equivalent widths 

For calibration clusters, the fixed-slope linear fits in the J^W 
vs. V — Vhb plane were performed for Y^W values on both 
systems (G09 or AD91 converted to G09), and the resulting 
(W) values were averaged. As described in Da Costa et al. 
(|2009p . the fits were computed only for stars with V— Vhb < 
0.2 because the (EW, V — Vhb) relation appears to notably 
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Table 3. Calibration Cluster Data 



ID 


RA (2000) Dec (2000) Rad Vel AV^ 


W"8542 e 


^8662 e 




(kms"""^) (mag) 


(A) 


(A) 



NGC3201 1 


760 


154.45472 


-46.41605 


486 


-1.02 


2.57 


0.14 


1.99 


0.11 


NGC3201 1 


2825 


154.43903 


-46.40734 


484 


-0.95 


2.44 


0.14 


1.86 


0.11 


NGC3201 1 


4251 


154.43742 


-46.40153 


487 


-2.12 


2.69 


0.12 


2.16 


0.10 


NGC3201 1 


5837 


154.43194 


-46.39489 


487 


-0.24 


2.19 


0.21 


1.65 


0.12 


NGC3201 1 


6876 


154.41135 


-46.39017 


501 


-1.60 


2.71 


0.14 


2.07 


0.08 



Notes. This table is available in its entirety in a machine-readable form in the online journal. A portion is shown here for guidance 
regarding its form and content. These data supersede those of Table 1 in Da Costa et al. (2009). 

Table 4. Program Cluster Data 









RA (2000) 


Dec (2000) 


Rad Vel 
(kms-i) 


(mag) 


W^8542 

(A) 


e 


W^8662 

(A) 


e 


Pyxis 


1 


1935 


136.97620 


-37.21546 


41 


-1.11 


2.31 


0.25 


2.15 


0.27 


Pyxis 


1 


3319 


136.97864 


-37.20316 


50 


-0.09 


2.52 


0.56 


2.08 


0.72 


Pyxis 


1 


3370 


137.02170 


-37.20279 


41 


-1.03 


2.78 


0.34 


1.67 


0.27 


Pyxis 


1 


4171 


137.04704 


-37.19520 


45 


0.02 










Pyxis 


2 


11646 


136.98251 


-37.27163 


67 


0.39 










Pyxis 


2 


13036 


136.97849 


-37.25888 


17 


-0.16 










Pyxis 


2 


14488 


137.01187 


-37.24609 


58 


-0.10 


2.27 


0.45 


1.36 


0.52 


Pyxis 


2 


16070 


136.99621 


-37.23422 


49 


0.00 


2.89 


0.47 


1.86 


0.61 



Notes. This table is available in its entirety in a machine-readable form in the online journal. A portion is shown here for guidance 
regarding its form and content. The M22 data in this table supersede those of Table 2 in Da Costa et al. (2009). 

flatten at lower luminosities, as predicted by models (see, 
e.g., Carrera et al. IMITl Starkenburg et al. [2UTD1) . This 
effect is shown particularly by the NGC6397 and MIO stars 
in Fig.[5]in which we show the calibration lines for the entire 
set of calibration cluster data. The YjW values adopted in 
the figure are the converted AD91 ones, but in both cases 
we found that the best-fitting slope is a = —0.627. 

This slope was then adopted to compute reduced equiva- 
lent widths for program clusters in both metallicity groups. 
In the case of NGC2808, Lynga 7 and Pal 7, which were 
measured with both methods, the same procedure as for 
the calibration clusters was adopted, and their (W') values 
are the average of the two results. 

Once reduced equivalent widths W' are computed, their 
average lyV') can be calibrated onto a metallicity scale us- 
ing standard clusters, as the next section explains. Finally, 
using the calibration relation, the metallicity of each star for 
all clusters in our database can be computed by converting 
its W' to [Fe/H], and a representative average metallicity 
can also be computed for each cluster. 

We note that, particularly in the differentially reddened 
clusters, the possibility that AGB stars are included in our 
"RGB" samples may lead to additional scatter in the (SW^, 
y — yBs) plane, and thus to increased uncertainty in the de- 
rived abundances. On the other hand this effect is expected 
to be small, as shown by Cole et al. (|2000p . 





1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
_ N6528 / R97 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
_ N6553 / R97 
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1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

7 N6528 / this paper 
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N6553 / this paper _ 
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5 
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^ ^--^^^ ♦ ^ 
^^^^ ♦ 

7 , , , , 1 , , , , 1 , , , , 1 , , , , 1 7 


.^^-^^^ ' 
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Figure 5. Comparison of YjW vs. V — Vhb using both our 
and R97 data. The R97 EWs were corrected using the fac- 
tor 1.117 found in Sect. 15.21 The solid lines are linear fits 
with fixed slope a — —0.627. When our data are plotted, the 
fitted lines to R97 data are also shown as dashed lines. A ra- 
dial velocity cut was imposed to R97a data to select cluster 
members: ARV> 38kms~i and ARV< 52.75 kms^^ define 
cluster members for NGC6528 and NGC6553, respectively. 
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4.3.1. NGC6553 and NGC6528 

All metallicity determinations in the literature put 
NGC6553 at a lower [Fe/H] than NGC6528, so it is some- 
what surprising that in this work we found NGC6553 
at a higher {W) than NGC6528. As Table E] indicates, 
(W') = 5.68 ±0.09 for NGC 6528 an d {W) = 5.84 ±0.03 for 
NGC6553. Rutledge et al. (Il997a|l found instead {W) = 
5.41 ± 0.14 and 5.13 ± 0.09 for the two clusters, respec- 
tively, which is consistent with their ranking in [Fe/H]. 
To understand the origin of this discrepancy we show in 
Fig. [5] the YjW vs. V — Vhb diagrams for the two clusters, 
using both R97 data and ours. Evidently, the problem is 
with NGC6528, for which we measure EWs that are ^0.6 
smaller than those of R97a. It should be recalled that mea- 
suring the strengths of CaT lines in high metallicity cluster 
stars is more difficult than at lower metallicities, because of 
the stronger presence of metal lines at lower temperatures. 
These contaminate line wings, which is where the increased 
strength resides, and they also make more difficult to de- 
fine the appropriate pseudo- continuum. The veiling caused 
by molecular species such as TiO and CN is also a po- 
tential concern at lower effective temperatures. The pres- 
ence of TiO depresses the pseudo-continuum and reduces 
the apparent strength of the CaT lines (see Olszewski et 
al. I199ip . and it might explain our small EWs. However, 
none of the NGC6528 members stars observed here or in 
Gullieuszik et al (2009) show any evidence for the presence 
of TiO bands in the part of the spectrum considered her^. 
Another possibility is that our V — Vhb is '^1 mag larger 
than that of R97a. There are no stars in common between 
the two datasets, however we have five stars in common 
with van den Bergh & Younger (1979; VY79), which tell 
us that we are fainter by 0.07 mag in the IZ-band. R97a 
used photometry from Ortolani et al. (|1992p . and found 
that it is fainter than that of VY79 by 0.05 mag in the 
same band. Therefore our zero-points should agree with 
those of R97a within 0.02 mag. In addition R97a adopted 
Vhb = 17.1 ± 0.1 from Ortolani et al. (1992), which is the 
same value that we used. Therefore it appears that a zero- 
point difference between us and R97a must be ruled out, 
but it is still possible that R97a measured six stars among 
which some had high extinction. For example excluding the 
faintest two stars from the R97a data, brings their average 
EW closer to our values by ^^0.2 A. Finally it is of course 
possible that either our Ca li line fits or those of R97a suf- 
fer from a still unidentified source of error. In this respect 
we can only remark that the trend defined by our data in 
Fig. [5] seems more 'natural' than that shown by the R97 
ones. Clearly an independent study is needed to solve this 
puzzle, possibly based on near infrared photometry like the 
works of Warren & Cole (pPUgj) and Lane et al. (PCTTTj) . 

4.4. Calibration of reduced equivalent widths 

The data needed to find the transformation between (W) 
and [Fe/H] values are summarized in Table [5] for each cali- 

^ In the case of NGC6553, five of the eighteen stars used 
to calculate the cluster mean velocity (see Table [5]) show def- 
inite indications of the presence of the A8440 AtIO bandhead. 
Consequently, these five stars have not been used in the line 
strength determinations for this cluster. The spectra of the re- 
maining thirteen NGC6553 stars show no sign of the presence 
of TiO bands. 



RMS = 0.125 




2 3 4 5 

<W'> G09 



Figures. ]Fe/H] on the C09 scale vs. (W) on the G09 
scale, from the data in Table [5] The solid curve shows the 
cubic fit to the data, and the ±RMS dispersion boundary 
is represented by the shaded area. Note that the error bars 
are smaller than the plot symbols. 

bration cluster we list its (W) and its error, and the metal- 
licity on three different scales. In Carretta & Gratton (jl997p 
and Kraft & Ivans ( 2003) there are no clusters more metal 
rich than [Fe/H] ~ —0.7 so our calibration relation is based 
on the Carretta et al. (|2009|) scale, nevertheless values on 
the other two scales are given for readers wishing to do com- 
parison studies in a more limited metallicity range. Using 
these data we found the following cubic calibration relation 
for the C09 scale: 

[Fe/H] = 0.0178 (W'f - 0.114 {W'f + 0.599 (W) - 3.113 

which has a RMS dispersion around the fit of 0.126 dex, and 
is defined in the [Fe/H] range from —2.33 dex to +0.07 dex. 
With this relation, W' values for each star in each pro- 
gram cluster can be converted into [Fe/H] values on the 
C09 scale, and average (W) can be converted into average 
metallicities, which are listed in TableEl The table also lists 
the uncertainty in our determinations which is the combi- 
nation of the statistical uncertainty in the (W) value, the 
calibration uncertainty (taken as the RMS dispersion about 
the fit) and, where necessary, inclusion of allowance for sig- 
nificant differential reddening that can affect the individual 
V — Vhb values. Also given in the Table are [Fe/H] values 
from HIO and their associated weight (higher numbers in- 
dicate more certain values) as well as the values and their 
uncertainties from Appendix 1 of C09. 

For the other two metallicity scales, the calibra- 
tion relations are [Fe/H]Kio3 = 0.496 x (W) - 3.369 
(RMS=0.063 dex) and [Fe/H]cG97 = 0.383 x {W) - 2.758 
(RMS=0.084 dex). The range of validity of these two re- 
lations is —2.45 dex to —0.67 dex for the KI03 scale, and 
-2.12 dex to -0.7 dex for the CG97 scale. 



4.5. Extrapolation of the calibration relation at low (W) 

When metallicity goes to zero, we expect that (W) goes to 
zero as well, while ]Fe/H] should go to — oo. Therefore our 
calibration relation cannot be extrapolated to arbitrarily 
low (W). This is illustrated by Fig. [71 where to our cali- 
bration clusters we have added two metal-poor stars from 
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Table 5. Input data for the (W) - [Fe/H] relation. 



NGC Alt. (W) e [Fe/H] e [Fe/H] e [Fe/H] src 

C09 CG97 KI03 



3201 




3.70 


0.03 


-1.51 


0.02 


-1.23 


0.05 


-1.46 


3 2 


4372 




2.33 


0.03 


-2.19 


0.08 






-2.29 


1 2 


4590 


M68 


1.85 


0.05 


-2.27 


0.04 


-1.99 


0.06 


-2.43 


1 2 


5904 


M5 


4.39 


0.02 


-1.33 


0.02 


-1.11 


0.03 


-1.32 


1 2 


5927 




5.29 


0.04 


-0.29 


0.07 






-0.67 


3 2 


6121 


M4 


4.36 


0.05 


-1.18 


0.02 


-1.19 


0.03 


-1.22 


3 2 


6171 


M107 


4.59 


0.05 


-1.03 


0.02 






-1.1 


1 2 


6254 


MIO 


3.37 


0.02 


-1.57 


0.02 


-1.41 


0.02 


-1.48 


3 2 


6397 




2.61 


0.04 


-1.99 


0.02 


-1.82 


0.04 


-2.12 


1 2 


6528 




5.68 


0.09 


0.07 


0.07 








1 3 


6553 




5.84 


0.03 


-0.16 


0.06 








3 2 


6752 




3.84 


0.02 


-1.55 


0.01 


-1.42 


0.02 


-1.46 


1 2 


6838 


M71 


5.09 


0.04 


-0.82 


0.02 


-0.7 


0.03 


-0.82 


3 2 


7078 


M15 


1.69 


0.04 


-2.33 


0.02 


-2.12 


0.01 


-2.45 


3 2 



Notes. The numbers in column src indicate the source of (W) measurements: (1) G09 paper, (2) AD91 method, (3) G09 method 
applied to clusters observed in this run. 



RMS = 0.122 




2 4 6 

<W> 



Figure 7. The metallicity - (W) relation is plotted here for 
calibration clusters (filled diamonds) plus two metal-poor 
stars from Starkenburg et al. (j2010l asterisks). Overplotted 
is a polynomial fit of degree n = 5, which shows that our 
calibration relation (dashed line) cannot be used outside 
the metallicity range defined by the standard clusters. 



Starkenburg et al. (j2010p . From their Fig. 9 we have es- 
timated the intercepts of their {T,{W), V — Vhb) relations 
at y — Vhb = 0, which give (W) for a range of metal- 
licitiefl In particular we have selected stars CD-38 245 
and HD88609, the two most metal-poor Galactic stars in 
Starkenburg et al. (PUTU)) . The two stars have [Fe/H]= -4.2 
and —2.9, so they can be used to see the trend of [Fe/H] vs. 
(W) beyond the range defined by our calibration clusters. 
Indeed Fig. [7] shows a downward trend for (W) — > 0, as 
expected. The position of the two lowest metallicity points 
is uncertain and their (W) are on a slightly different mea- 
surement system from that of G09, so they cannot be used 
to define a calibration relation valid over a larger [Fe/H] 



^ Note that for single stars an average (W) does not make 
sense, but for simplicity we use that symbol instead of W' . 



range. However for the sake of illustration the whole data 
set was fit with a polynomial of degree n ~ 5. The fig- 
ure makes clear that our fiducial relation cannot be used 
beyond the [Fe/H] range defined by our standard clusters. 



5. Discussion 

The average metallicities of the calibration and program 
clusters as derived from our CaT measurements are com- 
pared to those of HIO and C09 in the panels of Fig. [S] In 
both cases there is no indication of any systematic deviation 
with abundance and the mean offset is close to zero. The 
agreement indicates that all three measurements are on a 
consistent system, which is not unexpected. The scatter in 
the abundance differences about the offset seen in the fig- 
ure reflects both the uncertainty in our measurements and 
in the existing determinations, which generally come from 
a variety of heterogeneous sources. 

5.1. Individual Cluster Results 



In the following sections we discuss our new [Fe/H] val- 
ues for clusters that have issues with (a) data quality, (b) 
differential reddening or low statistics affecting the position 
of the HB, (c) abundance spread, (d) stars' membership, (e) 
field contamination, and (f) large [Fe/H] differences with 
the literature. A first assessment of possible abundance 
spreads was done by computing, for each cluster, the ratio p 
of the r.m.s. dispersion around the linear fit with fixed slope 
a = 0.627 A mag~^, to the average measurement error in 
T,W. This parameter is given in the last column of Table [51 
Clusters for which p > 1.5 are considered candidates for a 
metallicity dispersion, while clusters for which 1 < p < 1.5 
are considered marginal candidates. More details are given 
in the next sections. 
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Figures. Differences in metallicity between our determinations and those of the literature are plotted here against 
[Fe/H] values computed in this work. In the upper panels the comparison metaUicities are those of HIO, while in the 
lower panels they are those of Appendix 1 of C09. The dashed and dotted lines represent the average and ±cr of the 
metallicity differences, which were computed after retaining clusters with deviations smaller than Ixa from a preliminary 
fit (encircled single or double diamonds). The left panels show the 14 calibration clusters, and the right panels show the 
20 program clusters (17 for the comparison with C09/Appendix 1, see Table 16]). In the right panels single diamonds 
represent clusters with weight smaller than 3 in HIO. Error bars in the lower panels represent the quadratic sum of our 
errors and those of C09. 



5.1.1. Pyxis 

The spectra were obtained in the second night, which was 
partially cloudy, so the data are of much poorer quality 
than for the rest of sample. Three of the 8 stars that were 
used to determine radial velocity do not have line strength 
measures (the Gaussian fit failed). The given W' value is 
the weighted mean value, and the HB luminosity was deter- 
mined by considering stars within 2' of center, and taking 
the mean of red HB stars. 

Our value of —1.45 ± 0.14 dex is somewhat lower than 
that (—1.20) tabulated by HIO. It is, however, consistent 
with the spectroscopic determination of Palma et al. (|2000p 
who measured [Fe/H]zw = — 1.4 ± 0.1 from a spectrum 
at the CaT of a single Pyxis red giant. The CMD based 
photometric determinations of Irwin et al. (|1995[) . — 1.1 ± 
0.3, and Sarajedini & Geisler (|1996|) . — 1.2 ± 0.15 are also 
not inconsistent with this determination. With its dominant 
red horizontal branch morphology Pyxis is clearly a 'young 
halo' object (cf. Da Costa il995t Sarajedini & Geisler il996|l . 

5.1.2. Ruprecht 106 

At face value the dispersion in the HW values for the 9 
candidate members of this cluster (0.30 A) is significantly 
larger than the mean measurement errors (0.15 A), sug- 
gesting the possible presence of an intrinsic abundance 
range in this cluster. On further investigation, however, 
such a result seems unlikely. The large apparent disper- 
sion is driven by two stars: star 2_ 17214 which lies below 
the fitted -0.627 A mag"^ line and star 2_14496 which 
lies above it. If these two stars are excluded the remaining 
7 have dispersion in "EW of only 0.089 Aand a mean mea- 



surement error of 0.13 A. Star 2_17214 is the brightest in 
the observed sample and in the color-magnitude diagram 
(CMD) derived from our pre-imaging photometry it lies 
bluer and brighter than the red giant branch. This CMD 
location is confirmed when the stars observed spectroscopi- 
cally here are cross-identified with the photometric study of 
Sarajedini & Layden (OMT)) . Our star 2_17214 is SL37 and 
again it lies significantly to the blue of the well defined red 
giant branch in the Sarajedini & Layden color-magnitude 
diagram. It is hkely that 2_17214/SL37 is either an AGB 
or post- AGB star that is effectively "too bright" in V — Vhb 
for its YiW value: RGB stars with the same V — I color are 
~0.7 mag fainter and such a magnitude offset would place 
the star very close to the fitted line in the Y,W, V — Vhb 
plane. 

The second discrepant star, 2_14496 (SL519), is the 
faintest in our sample. With V — Vhb=0.18, it is close to 
the V — Vhb =+0.2 cutoff for inclusion in the fitting pro- 
cess, for fainter magnitudes the line strengths do not de- 
crease at the same rate as for more luminous stars. Using 
the Sarajedini & Layden (I1997P photometry this star has 
V — Vhb =0.25 and thus would be automatically excluded 
from the fit. Its line strength is not inconsistent with those 
of the slightly more luminous stars (cf. the fainter MIO 
and NGC6397 stars in Fig. |4]). We conclude there is no 
compelling reason to consider this star discrepant. Further, 
the other stars in our spectroscopic sample have locations 
in the Sarajedini & Layden (jl997p CMD consistent with 
our photometry and with locations on the RGB. Taken to- 
gether these arguments suggest that there is no evidence 
for any significant abundance spread in Ruprecht 106. This 
result is consistent with the earlier work of Da Costa, 
Armandroff & Norris (|1992p who observed 7 red giants in 
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Table 6. Metallicities of program clusters. 
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Notes. Column (1) lists the main and alternate cluster ID. Column (2) is the number of cluster members used in determining the 
{W) values given in Column (3). Column (4) gives the source of the {W) values: 1 for direct G09 measurements, 2 for transformed 
AD91 measurements and 3 for the average of both measurement techniques. Column (5) is the standard deviation of the (W) 
values. Column (6) is the [Fe/H] on the C09 scale which results from applying the abundance calibration. Column (7) is the total 
uncertainty in the [Fe/H] value: the rms sum of the abundance uncertainties from the standard error of the mean in {W}, from the 
calibration uncertainty, and from the effects of differential reddening where necessary. Columns (8) and (9) are the HIO abundance 
and weight while columns (10) and (11) are the abundance and error from Appendix 1 of C09 (clusters with a "1" in the Notes 
column only). Column (12) indicates if the effects of differential reddening are significant, details are given in the appropriate 
sub-section for these clusters. Column (13) gives the ratio of the r.m.s. dispersion around the linear fit with fixed slope, to the 
mean measurement error. In the last column we mark candidate clusters for metallicity dispersion with a solid arrow head, while 
marginal candidates are marked with an open arrow head. 



this cluster at the Ca ii triplelQ. They found an abundance 
of [Fe/H]= -1.69±0.05 dex fully consistent with that de- 
rived here. Moreover, the 7 stars observed by Da Costa et 
al. I|1992p have a dispersion in I]W about their adopted 
{EW, V - Vhb) relation of 0.16 Awith a mean measure- 
ment error of 0.25 A. Again this is consistent with a lack 
of any intrinsic abundance dispersion in this cluster. 

5.1.3. An abundance spread in NGC5824 



This object is a little studied, luminous {My — —8.85) 
globular cluster that lies 32 kpc from the Sun and 26 kpc 
from the Galactic center (HIO). After NGC 2419 for which 
My = -9.4 and Rq = 90 kpc, NGC5824 is the next most 
luminous distant globular cluster in the outer halo - M54, 
the nuclear star cluster of the Sgr dwarf galaxy has {My, 
Rg) of (-9.98, 18.8) while NGC 5024 has (-8.71, 18.4) us- 
ing the data tabulated by HIO. The only previous spectro- 
scopic abundance determination for this cluster is based on 

There is one star in common with our work: star 2205 of Da 
Costa et al. (1992) is star 16406 here. 



an integrated spectrum at the CaT obtained by ArmandrofF 
& Zinn (|1988p . This is the source of the abundances of 
[Fe/H]= -1.91 (weight 4) and [Fe/H]= -1.94 ±0.14 listed 
by HIO and C09, respectively. 

In the lower panel of Fig. [2] we show our T,W measures 
plotted against V — Vub for the 17 members observed in this 
cluster, together with the calibration lines for the standard 
clusters. The measurements were made on the AD91 sys- 
tem but have been transformed to the G09 system. The 
value of (W) for these stars leads to an abundance of 
[Fe/H]= -2.01 ± 0.13 for NGC5824, which is in excellent 
agreement with the earlier determination. However, while 
the effect is not as striking as it is for M22 (see Da Costa et 
al 2009 Fig. 4), it does appear in Fig. [5] that the dispersion 
of the NGC5824 T,W values is larger than would be ex- 
pected from the measurement errors alone. Specifically, the 
standard deviation about the fitted line is 0.31 A while the 
mean measurement error in "EW is 0.14 A, suggesting the 
presence of an internal abundance spread in NGC5824- The 
largest contribution to the measured dispersion in the T,W 
values comes from the stars 2 32429 and 1_8575, which 
have EW values 0.72 and 0.59 A., respectively, higher than 
that expected from the mean relation at their V — Vhb val- 
ues. These two stars are not distinguished from the other 
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Figure 9. Data for stars observed in NGC5824. The up- 
per panel shows the observed hehocentric velocity plotted 
against V—Vhb- The dotted line is the cluster velocity given 
by the Harris catalogue, the dashed line is the mean veloc- 
ity for the 17 stars selected as cluster members. The sole 
star identified as a radial velocity non-member is plotted as 
an x-sign. The lower panel shows 'EW against V — Vhb for 
the same stars. The dot-dash line is the fit of the calibra- 
tion relation to the cluster member points (filled circles). 
Shown also are the calibration lines for the standard clus- 
ters M15, NGC4372, NGC6397, MIO, M5 and M71 (in order 
of increasing T,W) . The significant scatter in the NGC5824 
member data suggests the presence of an intrinsic metallic- 
ity dispersion. 



NGC5824 stars observed in terms of radial velocity (see 
upper panel of Fig. distance from the cluster center, 
or location in the color-magnitude diagram. Further, the 
T,W measurement error for star 2_32429, 0.17 A, IS con- 
sistent with the mean for the other stars. The value for 
1_8575 is somewhat higher at 0.23 A, but that is not sur- 
prising given it is the faintest star in the observed sample. 
If this star is excluded the dispersion in the IIW values is 
reduced only marg inally from 0.31 Ato 0.28 A, while the 
mean error in Y^W is unchanged at 0.14 A. We conclude 
that despite the small number of stars observed there is 
strong evidence for intrinsic line strength variations in our 
NGC5824 sample. For comparison, the equivalent numbers 
for M22 are (0.30, 0.11) with a sample of 41 stars, while for 
M54, a cluster also known to have an internal abundance 
spread, e.g. Carretta et al. (|2010l hereafter CarlO), they are 
(0.25, 0.16) for a sample of 15 stars (see below). These are 
the only three clusters in our entire sample where there is 



substantial evidence for the presence of an intrinsic [Fe/H] 
dispersion. Thus NGC5824 now joins the small number of 
globular clusters where such intrinsic [Fe/H] dispersions are 
known 

If we subtract the contribution from the equivalent 
width measurement error and apply the abundance cali- 
bration, then the corresponding intrinsic abundance dis- 
persion in NGC5824 is crint([Fe/H]) = 0.12 dex, which is 
comparable to our results for M22 ((Tint ([Fe/H]) = 0.15 
dex, Da Costa et al. 2009). A dispersion of this order 
places NGC5824 at a consistent location for its luminos- 
ity in Fig. 7 of CarlO. With a sample of only 17 stars 
the form of the abundance distribution is not well con- 
strained but the impression from Fig. [HI is that it may 
be similar to those of M22 (Da Costa et al., 2009) and 
wCen (e.g., Johnson & Pilachowski 12010]) with a steep rise 
on the metal-poor side to a peak at ]Fe/H[= —2.06 and 
a broader tail to [Fe/H[Ri —1.7 dex. The median abun- 
dance is [Fe/H]= —2.02 dex. We have accepted programs at 
Gemini-South with CMOS and at the VLT with F0RS2 to 
substantially increase the number of NGC5824 red giants 
with abundance determinations. 

The discovery of a probable intrinsic abundance spread 
in NGC 5824 is particularly intriguing as Newberg et al. 
(12009 *) have recently suggested that this cluster is possi- 
bly the former nucleus of a dwarf galaxy whose tidal dis- 
ruption is responsible for the stellar stream known as the 
Cetus Polar Stream (Newberg et al BUOg)) . The Cetus Polar 
Stream is a low metallicity ([Fe/H[« —2.1) tidal stream ap- 
proximately 34 kpc from the Sun. A connection between the 
stellar stream and NGC5824 would strengthen the hypoth- 
esis that there exists a set of globular cluster-like stellar sys- 
tems, characterized by the presence of internal [Fe/H] abun- 
dance ranges and higher than average luminosity, which 
have their origin as former dwarf galaxy nuclei (e.g., wCen) 
or dwarf galaxy central star clusters (e.g., M54) and which 
are distinct from "regular" globular clusters (e.g. Da Costa 
et al. 2009, CarlO). 

5.1.4. Lynga 7 

Stars within a radius of 1' from the cluster center were 
initially selected from the photometry dataset to define 
the Vhb value and to define the cluster locus in the 
(SW, V — Vhb) plane. This process and the individual ra- 
dial velocities indicated 5 stars as probable members. A 
further 4 stars were then added to the sample based on ra- 
dial velocities and line strengths consistent with the prob- 
able members. However, in calculating the value of (W), 
the brightest star (1_4175) was excluded as its spectrum 
clearly reveals the presence of TiO absorption bands that 
affect the strength of the CaT features. We note also that 
our CMD clearly shows the effects of differential redden- 
ing across the cluster in that the blue end of the horizontal 
branch is ^ 0.5 mag brighter than the red end. We used 
the midpoint value for Vhb and note that a ±0.13 mag un- 
certainty in Vhb (i one quarter of the full range) gives an 
additional uncertainty of ±0.06 dex in the abundance deter- 
mination. Our determination of [Fe/H] C09 = — 0.57 ±0.15is 
larger than that tabulated by HIO but is entirely consistent 



From a preliminary analysis of high resolution spectra of 
three stars from our sample Villanova & Geisler (in preparation) 
find a significant [Fe/H| difference, much larger than the errors. 
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with the previous spectroscopic determination of Tavarez & 
Friel ([TM5)) who found [Fe/H] = -0.62 ±0.15 from an anal- 
ysis of Fe line strengths in moderate resolution blue spectra 
of 4 Lynga 7 red giants. The relatively high abundance is 
consistent with the interpretation of Lynga 7 as a (thick) 
disk globular cluster (cf. Ortolani et al. I1993p . 

5.1.5. NGC6139 

Our [Fe/H] value is the first spectroscopic determination 
for this cluster, previous determinations derive from the 
integrated light Qag photometry index of Zinn (|1980p . or 
from the colour of the red giant branch in the CMD. 

5.1.6. Terzan 3 

Three probable members included in the calculation of the 
cluster radial velocity were excluded from the W' determi- 
nation as they are clearly fainter than V — Vhb = -f 0.2. 
As for Lynga 7, this cluster suffers from notable differential 
reddening with the blue end of the predominantly red hor- 
izontal branch being approximately 0.6 mag brighter than 
the red end. The midpoint was used as the V — Vhb value. 
Again assuming the uncertainty in this value is ± one quar- 
ter of the full range, the corresponding uncertainty in the 
derived abundance from the effects of differential redden- 
ing is 0.06 dex. The only previous spectroscopic determina- 
tion of the abundance of Terzan 3 is that of Cote (|1999p . 
who found [Fe/H]= —0.75 ± 0.25 from an analysis of the 
strengths of strong Fe 1 lines on low S/N high resolution 
spectra, which is 0.33 dex higher than our value. Given the 
lack of well determined abundances for the clusters used 
to calibrate his line strength - abundance relation, the dis- 
crepancy with our determination is not significant. 

5.1.7. NGC6325 

This cluster is marked in Table IH] as having a marginal 
metallicity spread. The dispersion of T,W around the fit vs. 
V — Vhb is inflated by one star only, 2_62566; which is also 
the faintest in the sample. If that star is removed from the 
fit, then the ratio of the rms dispersion around the fit to the 
average error in EVF becomes p — 0.92. We conclude that 
there is no evidence of a metallicity dispersion in NGC6325. 

5.1.8. HPl 

This cluster is projected against a dense Galactic bulge 
field and consequently isolating candidate members is not 
straightforward. We began by considering the CMD for 
stars within 23" of the cluster center (cf. Ortolani et al. 
I1997P to identify the blue horizontal branch of the cluster 
and thus establish Vhb- The resulting (SW, V — Vhb) and 
radial velocity data then allowed the selection of 6 proba- 
ble members within 1' of the cluster center and a further 2 
probable members at larger radial distances. Of the 6 more 
central stars, 5 are confirmed as members also by the proper 
motion analysis of Ortolani et al. (|201ip . while the other 3 
are outside their explored area. The abundance determina- 
tion is based on all 8 candidate members but is unchanged 
if only the six inner stars are considered. 

Our abundance [Fe/H]= — 1.32±0.14 is somewhat lower 
than the value -1.0 ± 0.2 given by Barbuy et al. (|2006p 



which was based on high dispersion spectra of two stars. 
Since [Ca/Fe]=0.03 was obtained in Barbuy et al. ((2006)) . 
the lower [Fe/H] value found here could not be explained 
by a Ca overabundance, as it may be the case for other 
clusters. The stars analyzed by Barbuy et al. (|2006p . HPl- 
2 and HPl-3, correspond to our stars 1_6931 and 1_4996 
and these are among our probable cluster members. Our 
lower abundance alleviates to some extent the disparity be- 
tween the blue horizontal branch morphology of this cluster 
and the relatively high abundance found by Barbuy et al. 
(pnU51) . We note also that recently Valenti et al PUTP)) de- 
termined an abundance of [Fe/H]cG97 = —1.12 ± 0.2 dex 
for HPl based on the infrared colors of the red giant branch 
stars. 

In Table [S] the value of p derived from the 8 candi- 
date members suggests the possible existence of an intrin- 
sic abundance spread in this cluster. However, we do not 
think this is likely to be the case. Rather we suggest that a 
combination of differential reddening, increased photomet- 
ric errors due to the crowded nature of the field, and the 
possible inclusion of AGB stars in the candidate member 
sample, have all contributed to the scatter about the fitted 
line marginally exceeding that expected from the measure- 
ment errors. Certainly the CMD for stars within 1' of the 
cluster center shows considerable scatter, and without a 
larger sample of confirmed members it is difficult to con- 
vincingly identify the cluster RGB and possible AGB se- 
quences. Spectroscopic observations of a larger sample of 
cluster members would of course also lead to stronger con- 
straints on the presence or absence of an intrinsic abun- 
dance spread in this cluster. 

5.1.9. NGC6380 

The brightest cluster member (1_3509) was not included 
in the calculation of {W) as the spectrum shows clear signs 
of TiO absorption affecting the location of the pseudo- 
continuum in the vicinity of the CaT lines. This cluster 
is also affected by significant differential reddening with 
the V magnitude difference between the blue and red ends 
of the predominantly red horizontal branch differing by 
'-^ 0.7 mag. We again adopted the midpoint for Vhb- As 
for Lynga 7 and Terzan 3 above, if we take ±0.18 mag 
as the uncertainty in the V^b value (i.e. ± one quarter of 
the full range), then there is an additional uncertainty of 
±0.09 dex in our abundance determination. 

5.1.10. NGC6440 

This cluster is another that suffers from notable differential 
reddening. Using the same approach as for the other clus- 
ters with differential reddening, the estimated uncertainty 
in the adopted Vhb magnitude is ±0.14 mag which leads to 
an additional uncertainty in our [Fe/H[ value of ±0.07 dex. 

5.1.11. NGC6441 

As for NGC6440 this cluster also suffers from differential 
reddening. Using the same approach as for the other differ- 
entially reddened clusters the uncertainty in the adopted 
Vhb value is ±0.15 mag. This leads to an additional uncer- 
tainty of ±0.07 dex in the [Fe/H[ value which is included 
in the total uncertainty listed in Table |51 
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5.1.12. NGC6558 

The brightest cluster member (1_4901) in our sample 
was excluded from the calculation of (W) as the spec- 
trum shows clear TiO absorption. Our abundance, [Fe/H]= 
— 1.03 ± 0.14, is in excellent agreement with that, —0.97 ± 
0.15, found by Barbuy et al. (|2007|l from an analysis 
of VLT/FLAMES spectra of 5 red giant members. Since 
[Ca/Fe]=+0.05 was derived, in this case the metallicity 
[Fe/H]w[Ca/H]. In contrast, HIO and C09 hst a lower abun- 
dance for this cluster. The origin of those values lies with 
the integrated photometry and integrated spectroscopy 
Q39 values discussed in Zinn & West (jl984|) . Given our 
more reliable determination and the agreement with the 
results of Barbuy et al. (|2007p our abundance is preferable. 
The presence of TiO in the spectrum of the brightest red gi- 
ant in our sample is also indicative of a relatively high abun- 
dance. We note that our stars 1_3494, 1_4901 and 1_5672 
correspond to Barbuy et al. (2007) stars F97, F42 and B73. 
With this relatively high abundance and blue horizontal 
branch morp hology (Rich et al. HMEl Barbuy et al. [WTI 
Barbuy et al. I2009P the cluster is similar to HPl. As dis- 
cussed in Barbuy et al. (j2007|) . these moderately metal-rich 
globular clusters with blue horizontal branch morphologies 
may have been amongst the first objects to form in the 
Galactic bulge. 

5.1.13. Pal7 

This cluster is also among those affected by differential red- 
dening. We have taken the uncertainty in the Vhb value 
for the dominant red horizontal branch in this cluster as 
±0.16 mag, which leads to an additional uncertainty of 
±0.08 dex in our abundance determination. The only other 
abundance estimate for this cluster is that of Cote (J1999') 
who derived an abundance of [Fe/H]= —0.75 ± 0.25 from a 
similar analysis process to that described in the results for 
Terzan 3. As noted in the discussion for that cluster, given 
the lack of well determined abundances for the clusters used 
by Cote (jl999|) to calibrate his line strength relation, we 
consider our determination significantly more reliable. 

5.1.14. NGC6569 

Like many of the clusters in our sample, NGC6569 is also 
affected by differential reddening. The uncertainty in the 
Vhb value of this red horizontal branch cluster is ±0.20 mag 
which leads to an additional uncertainty in our abundance 
determination of ±0.09 dex. Our determination is some- 
what less than the values listed by HIO and C09. These 
have their origin in the Q39 value given by Zinn & West 
(|1984p . which was derived from integrated photometry and 
integrated spectroscopy of the cluster. Recently, Valenti et 
al. (2011) Hst [Fe/H]= -0.79 ± 0.02 (internal error only) 
for this cluster based on high resolution near-IR spectra of 
6 cluster red giants. This value is only slightly more than 
1 sigma different from our determination. There are no stars 
in common between our sample and that of Valenti et al. 

HMD). 

5.1.15. M22 

Our results for this cluster have been discussed in Da Costa 
et al. (j2009p . Here we note only that a different mean abun- 



NGC 6715 (M54) 




1 -1 -2 -3 

Figure 10. Y^W against y — Vhb for M54 stars compared to 
calibration clusters, showing a clear metallicity dispersion. 
Five-point asterisks identify the three metal-rich stars that 
are likely members of the Sagittarius dwarf galaxy popula- 
tion. 

dance results from the use of the C09 abundance scale 
rather than that of KI03. 



5.1.16. M54 

M54 is a luminous globular cluster superposed on the nu- 
cleus of the Sgr dwarf galaxy (e.g. Bellazzini et al. 2008]). 
The most recent detailed study of this system is that of 
Carlo who analyzed high dispersion VLT/FLAMES spec- 
tra for somewhat more than 100 stars in the central field 
of Sgr. They associate the bulk of their stars with M54 
finding a mean abundance of —1.56 ± 0.0^ and an intrin- 
sic abundance dispersion of 0.19 dex (76 stars). The re- 
mainder of their stars are relatively metal- rich ([Fe/H]> 
— 1.2, see Fig. 4 of CarlO) and CarlO associate them with 
the Sgr nucleus population finding a mean abundance of 
— 0.62±0.07 dex with a large intrinsic dispersion of 0.35 dex 
(27 stars). 

In Fig.[TU]we plot our T,W measures against V — Vhb for 
the stars observed in our M54 field. The measurements were 
made on the AD91 system but have been transformed to 
the G09 system. These data show 3 stars (1_5389, 2_20946 
and 2_22704) with relatively large values of T,W while the 
remaining 15 have weaker line strengths. We suggest that 
these 3 stars, which have abundances of [Fe/H]= —1.24, 
—0.97 and —1.21, respectively, belong to the Sgr population 
rather than M54. We note that it is not surprising that 
we did not find any more metal-rich members of the Sgr 
population as the original target selection was deliberately 
biased towards the M54 RGB population (see CarlO Fig. 3). 
The remaining 15 stars, assumed to be members of M54, 
have a mean abundance of —1.73 ± 0.13 dex. This is in 
reasonable agreement with the mean abundance —1.56 ± 
0.02 given in CarlO. 

As was found for M22 (Da Costa et al. 2009), and for 
NGC5824 above, the dispersion in the T,W values in Fig.fTUl 



^ Note that Carretta et al. (ppng]) listed [Fe/H]= -1.44 ±0.07 
for M54. 
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is notably larger than that expected on the basis of the 
uncertainties in the Y,W values. The dispersion about the 
fitted line in Fig.[Tn]is 0.25 A while the mean uncertainty in 
the T.W values is 0.16 A. Subtracting this error contribu- 
tion in quadrature then gives an intrinsic dispersion in SVF 
of 0.19 A, which translates to an intrinsic abundance dis- 
persion (TiNT([Fe/H]) of 0.09 dex. This is smaller than the 
0.19 dex intrinsic dispersion found by CarlO. We have no 
straightforward explanation for this discrepancy other than 
to note that the CarlO sample is five times larger and may 
therefore more fully probe the extremes of the underlying 
distribution. Uncertainties in the methods used to fix T^s 
might explain shifts in the absolute metallicities, but they 
should not affect relative abundances, and thus abundance 
dispersions. 

There are only 3 stars in common between our study 
and that of CarlO. Our M54 stars 1_1982, 2_10700 
and 2_ 17683 correspond to CarlO M54 stars 38004437, 
38004707 and 38009987, respectively, while none of our 
three postulated Sgr stars are in the CarlO list. For the 
three stars in common, the differences in the [Fe/H] val- 
ues, in the sense (our value - CarlO) are 0.15, —0.21 and 
—0.08 dex for a mean difference of —0.05 dex with a sigma 
of 0.18 dex. The mean difference is in the same sense as the 
difference in the mean abundances (—0.10 dex). The dis- 
persion, however, is somewhat larger than expected given 
that our individual abundance determinations have an un- 
certainty of order 0.11 dex and those of CarlO nominally of 
order 0.02 dex. Without a larger sample of common objects 
this question cannot be investigated further. 

5.1.17. Terzan 7 

Five probable cluster members were not included in the cal- 
culation of (W) as they are fainter than V — Vhb= +0.2. 
Their line strengths are, however, consistent with those of 
the more luminous members. Further, two additional prob- 
able members were excluded from the {W) calculation as 
their line strength measures were particularly uncertain. 

5.1.18. NGC7006 

Two probable cluster members were excluded from the (W) 
calculation as they are fainter than V — Vhb = +0.2. 
Their line strengths are nevertheless consistent with those 
of the more luminous members. For the remaining 18 
stars the dispersion in 'SW about the fitted line of slope 
—0.627 A mag^^ is 0.26 while the mean measurement error 
is 0.21 A. While this might be considered as evidence for the 
presence of a small intrinsic metallicity spread, we do not 
believe this to be the case. The increased dispersion in Y.W 
is driven by 3 stars that lie below the fitted line. Inspection 
of the location of these nominally discrepant stars in the 
CMD derived from our pre-imaging photometry indicates 
that all three stars are likely AGB stars. At the same color 
AGB stars are brighter in V — Vhb than RGB stars and 
this magnitude offset can reach 0.5 mag or more for stars 
on the lower part of the AGB, as is the case here. It results 
in AGB stars lying to the right of RGB stars of similar 
color in the (SM^, V — Vhb) plane mimicking a lower abun- 
dance. If these three AGB stars are excluded from the fit, 
the dispersion is reduced to 0.20 A and equals the mean 



RMS = 0.093 




2 4 6 

W'(R97) 

Figure 11. In the upper panel reduced equivalent widths of 
calibration clusters on our scale are plotted against those 
of R97. The data were fit with a linear relation passing 
through (0,0), and a la rejection was also applied after 
a preliminary fit. This leaves the clusters represented by 
encircled dots. A free linear fit to the whole sample is rep- 
resented by the dotted line. In the lower panel the difference 
of our (W) minus the converted R97 ones is shown, and 
the average difference ±cr are shown by the solid and dotted 
lines, respectively. 



measurement error. We conclude that there is no evidence 
in our data for an intrinsic abundance spread in this cluster. 

5.2. A merged catalog of globular cluster metallicities 

To enlarge the sample of clusters with homogeneous metal- 
licities, we took the reduced equivalent widths from R97 
and converted their (W) values to our (W) scale. As 
Fig. [m shows, R97 widths can be converted to our scale 
by the simple relation (W^')qo9 — cl x (W^')r97j where 
a — 1.117+0.043, which was obtained by comparing {W) of 
calibration clusters. For these clusters, after a la rejection, 
the average difference + sigma between our (W) and the 
converted R97 ones is 0.0003 + 0.0984 dex (lower panel of 
Fig. [TT|) . Apart from calibration clusters, we have four ob- 
jects in common with R97 (NGC2808, Rupl06, NGC6715, 
and Terzan 7) which can be used to test the transforma- 
tion between the two systems. After excluding NGC6715 
because of the metallicity dispersion, the comparison shows 
that for the remaining three clusters our (W) are on aver- 
age 0.15 ± 0.18 A larger than the converted R97 ones. This 
offset might seem large, but on the other hand there are 
also two calibration clusters with a similar deviation (see 
lower panel of Fig. [TT|l . We therefore expect that a larger 
set of common clusters would yield a null average differ- 



16 



Saviane et al.: Homogeneous Metallicities and Radial Velocities for Galactic Globular Clusters 



0.5 



> -0.5 



0.5 



-0.5 



<A[Fe/H]> = -0.01 ±0.06 




<A[Fe/H]> = -0.01 ±0.05 




-2 -1 

[Fe/H], 



Figure 12. Same as Fig. [5] for R97 clusters. Note that the 
[Fe/H] for the most metal-rich data point was extrapolated 
beyond the range of our calibration. 



ence between the G09 and the R97 converted (W) as for 
the calibration clusters. 

Once the (W) of R97 were transformed onto our sys- 
tem, we could apply the same calibration used to transform 
our {W) values to metallicity values on the C09 scale. The 
new R97 metallicities together with [Fe/H] from HIO, and 
from Appendix 1 of C09, are given in Table [71 Figure [T2l 
shows a comparison of our new [Fe/H] values to those of 
HIO and Appendix 1 of C09. The same comparison for our 
program clusters was shown in Fig. [8l and an inspection of 
the two figures reveals that the metallicity differences in the 
case of R97 clusters show a scatter which is about half that 
for our program clusters. This is not unexpected since the 
R97 compilation of equivalent widths is at the base of all 
recent abundances scales, so for those clusters the internal 
scatter is small. In the future the twenty clusters studied 
here will enter the set of objects with well-determined and 
homogeneous metallicities. 

5.3. Impact of the new abundances on the system of GGCs 

A common way to classify GGCs is to use the [Fe/H[ vs. 
HB type diagram, first introduced by Zinn (jl993l Z93). 
Such diagram for our clusters is shown in Fig. 1131 where 
the HB type is taken from Mackey & van den Bergh (2005). 
The figure shows that in 50% of the cases our new metal- 
licity values do not change the clusters' classification sig- 
nificantly, but for the other 50% of objects the change in 
[Fe/H[ implies a change in age of ^ 1 Gyr. In particu- 
lar NGC7006, NGC6569 and NGC6715 move to younger 
ages, and NGC6558 moves to older ages, if our new abun- 
dances are adopted. The cases of HPl and NGC6380 are 
instead more ambiguous, because their change in age de- 
pends on which comparison abundances are used. With our 
new metallicity determination the location of HPl in the 
Z93 diagram moves in the direction of younger ages when 
compared with its location using the HIO metallicity es- 



timate, and in the direction of older ages when compared 
to its location with the C09/ Appendix 1 abundance esti- 
mate. Our new metallicity for NGC6380 is similar to that 
of C09/ Appendix 1 implying a null age change, but the clus- 
ter becomes older if our new [Fe /H] is used instead of that 
of HIO. There are three clusters that do not have an entry 
in C09/ Appendix 1, so based on the change with respect 
to the HIO metallicity only, our new ]Fe/H] values would 
mean that both Pyxis and Terzan 3 become younger, while 
Lynga 7 becomes older. 

In the classification of Z93 the division between so-called 
"young halo" and "old halo" objects starts at an isochrone 
intermediate between the oldest two of Fig. [T3| (see, e.g. 
Fig. 5 of Mackey and Gilmore l2004p . so the greatest impact 
of our new abundances is for NGC 6715 (M54) which moves 
into the "young halo" area of the diagram. For clusters 
with HIO metallicity only, Lynga 7 becomes "old halo" and 
Terzan 3 becomes "young halo". Pyxis was already young 
halo, but it becomes the youngest cluster in our sample, 
together with Rupl06. 

Very recently Dotter et al. (j201ip have provided age es- 
timations for six GGCs at Galactocentric distances i?GC > 
15 kpc, including three clusters of our sample. They found 
that Rupl06, and Pyxis are 1-2 Gyr younger than inner 
halo GCs with similar metallicities, and that NGC7006 
is marginally younger than its inner halo counterparts. 
Figure[T3|is consistent with Dotter et al. in terms of relative 
ages, but the ages predicted by Rey et al. (j2001[) isochrones 
are lower by 1 Gyr. 

6. Conclusions and outlook 

The main result of this work is presented in Table [HI which 
gives the average metallicity on the C09 scale for our twenty 
program clusters, most of which are globular clusters in the 
distant halo or the Galactic bulge, with no or scanty spec- 
troscopic observations. Together with 68 clusters in Table [7| 
(72—4 in common), they add up to 88 clusters having metal- 
licities on the same scale, which is currently the largest set 
of objects with [Fe/H[ measured with the same method. In 
addition Table [2] gives also homogeneous radial velocities 
for our program clusters, which can be important for kine- 
matical studies of the GGC system. 

A comparison with literature metallicities revealed that 
our clusters had [Fe/H] values with uncertainties at the 
level of 0.2 dex RMS, because of inconsistent zero-points 
among different authors (Fig. [8]) . In particular the largest 
revisions of [Fe/H[ values happen for NGC6569, NGC6715, 
NGC6558, HPl, NGC6380, and NGC7006. When analyzed 
in the light of the [Fe/H[ vs. HB-type diagram, most 
clusters appear younger with the new metallicities. For 
NGC 6715 (M54) the younger age means that it now be- 
longs to the "young halo" area of the diagram. 

For all clusters we measured metallicities for an average 
of 10 stars, which let us search for metallicity dispersions. 
The M22 dispersion was discussed in Da Costa et al. (2009), 
and we confirmed the [Fe/H] dispersion of M54 with 15 
stars. In addition with 17 stars measured, we detected a 
probable dispersion also in NGC 5824. 

Some of our targets were included as 'special' objects, 
in particular Terzan 7 as Sgr cluster, NGC7006 as sec- 
ond parameter object, and NGC6325, NGC6356, NGC6440, 
NGC6441, and HPl as bulge members. The new Terzan 7 
metallicity is not so different from that in HIO and C09, 
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Figure 13. The metallicity of our program clusters on the C09 scale is plotted here against HB-type from Mackey and 
van den Bergh (|2005p . Isochrones are from Rey et al. (|200ip . and are separated by 1.1 Gyr from top to bottom. The 
oldest isochrone gives the age of clusters at i? < 8 kpc (Rey et al. 2001). Pal 7 is not present in Mackey and van den Bergh 
l|2005p . so it is not plotted here. The arrows connect the position of the cluster if another [Fe/H] estimate is used, to the 
position determined by our metallicity scale. The other scale is HIO on the upper panel, and Appendix 1 of C09 in the 
lower panel. Consistent changes for both scales are those of NGC7006, NGC6569, and NGC6715 which become younger, 
and that of NGC6558 which becomes older. Pyxis, Terzan 3, and Lynga 7 appear only in the upper panel because they 
have no [Fe/H] estimates in C09. Comparing our abundances to HIO only, Lynga 7 becomes older, while the other two 
become younger. In fact Pyxis would now be as young as Rupl06. 



so we confirm that it follows the age-metallicity relation 
of the Sgr dwarf. NGC 7006 was the first cluster where 
a violation of the HB-morphology vs. [Fe/H] relation was 
found, by Sandage & Wildey (1967). Its HB is too red for 
its metallicity. Our new metallicity is lower than literature 
values, so the violation becomes even more severe, although 
in the framework of Z93 it would be explained by an age 
lower than fiducial "old halo" clusters, as shown by Fig. [T31 
For bulge clusters our study, besides providing in many 
cases the first spectroscopic metallicity measurements of in- 



dividual stars, allows a useful pre-selection of cluster mem- 
bers, to be followed up with high-resolution studies. By 
comparing Tables [1] to [5] we can see that only 20% up to 
58% of stars observed in these clusters were confirmed as 
members, so our work will indeed allow future campaigns to 
save substantial observing time, not having to worry about 
foreground/background contamination. 

Finally the detection of a possible metallicity dispersion 
in NGC5824 deserves confirmation with a larger database 
of stellar [Fe/H], which we are in the process of collecting 
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both at the VLT and Gemini. Interestingly, clusters where 
dispersions have been detected are all possible nuclear clus- 
ters of dwarf galaxies in different stages of merging with 
the Milky Way, and they also have multiple stellar popula- 
tions (see, e.g., Piotto [^009,1 . In this context, with our new 
metallicity M54 now falls on the same low-age isochrone of 
NGC5824 in the Z93 diagram (see Fig.[T31), so this feature 
might become another distinctive trait of at least some nu- 
clear clusters. In this respect it is interesting to note the 
work of Paudel et al. (|201ip : they determined the ages of 
the stellar populations of 26 early-type dwarf galaxies in 
the Virgo cluster, and they found that in most cases (21 
galaxies) their nuclei are younger than the main bodies. 

We expect that our database of homogeneous metallic- 
ities for ~ 90 clusters (and radial velocities for 20 clus- 
ters) will be the foundation of many future studies of age- 
metallicity and age- kinematics relations. We plan to collect 
data on the 29 remaining clusters to complete our project. 
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Table 7. Metallicities for globular clusters from R97 on the Carretta et al. (2009) scale 
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Figure A.l. The sum of equivalent widths for the strongest 
two CaT Hues measured with Gaussian-only fits are 
plotted here against the same quantity measured with 
Gaussian+Lorentzian fits. The latter are taken from the 
G09 paper, for stars in the clusters M5 (+ signs), NGC4372 
(filled triangles), NGC6752 (filled 5-point stars), NGC6397 
(open 5-point stars), NGC 6171 (open squares) and NGC 
4590 (6-point stars). The interpolating line has slope / — 
1.147, and it was obtained by minimizing residuals and im- 
posing the passage through the origin. 

Appendix A: Transformation between the AD91 
and G09 measurements 

We define the reduced equivalent width for a single star as 
W ^ S - a{y -VYSi) where S = W^8542 + WgeeH So we 
can write for the (G)09 and (A)D91 systems: 

W^ = Sg- ao {V - 1/hb) 

W^^Sa- OA {V - Fhb) 
And the averaged reduced equivalent widths are: 

(Wi^) = [S5g - ac S(V - Vhb)]/N 

{W'a) = [S^A - OA ^{V - Vhb)]/N 

By comparing measurements of individual S made with the 
G09 or the AD91 method, we found that Sq = f Sa, and by 
imposing a linear fit that is passing through (0, 0) we found 
for clusters of the Leo I paper / = 1.147 (see Fig. IA.l| ) and 
for clusters in this run / = 1.090 (see Fig. IA.2| l. So then in 
the G09 system: 

(Wl,) = [f E5a - ac - VnB)]/N 
Then we have 

S^A = N (W^) + aA E(y - Vhb) 

^'^ We use S here instead of TlW because in this way the for- 
mulas are easier to read. 
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Figure A. 2. Same as Fig. lA.ll for the stars observed for 
this work. The symbols identify NGC3201 (open circles), 
MIO (open triangles), M4 (open diamonds), M15 (x signs) 
and M71 (open squares). The fitting line has slope / = 
1.090 and passes through the origin. 

And with a substitution 

(W^g) = [/ (WI) + "A She) - ac ^kb]/N 

where for simplicity we wrote She = S(V" — Vhb)- 
Rearranging and simplifying a bit more: 

N (Wi^) = fN (W^) + (/ aA - aG)SHE 

Which gives {Wq) as a function of / (with / — 1.090 or / = 
1.147), the number of stars N measured for each cluster, 
the averaged reduced equivalent width in the AD91 system, 
the slopes a a and oq of the S' vs. V — Vhe relations in 
the two systems, and the sum of the magnitude differences 
with respect to the HB, She- The difference in slopes / is 
probably due to the fact that older spectra had a better 
resolution (slit of O'.'S instead of 1"), so it might be that 
the Gauss-only fit does a worse job with line wings in that 
case. 
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